Field angle-dependent thermal conductivity in nodal superconductors 
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We apply a semi-classical method to the problem of field angle-dependent oscillations of the 
density of states and thermal conductivity for nodal superconductors and apply our results to the 
superconductor PrOs4Sbi2. The oscillatory contributions to the thermal conductivity for all possible 
point node configurations for a superconductor with Th symmetry are calculated. It is found that 
experimental results are best accounted for by nodes in the off-axis directions [isin^OjO, ± cos^o], 
which are associated with the time- reversal breaking, triplet paired phase with symmetry D2{E). 

PACS numbers: 74.20.-z, 71.27.+a, 71.10.-w 



I. INTRODUCTION 

The low temperature thermodynamic properties of 
unconventional superconductors are governed by nodal 
quasiparticles. In the presence of a small magnetic field 
Hci < -ff <C Hc2, it was shown by Volovik^ that the 
dominant contribution to the density of states (DOS) for 
superconductors with line nodes comes from delocalized 
quasiparticles, in contrast to s-wave superconductors in 
which the DOS is dominated by quasiparticles localised 
inside vortex coresJ^ Volovik argued that the delocalised 
states experience a semi-classical adjustment to their en- 
ergy due to the magnetic field which is expressed as a 
Doppler shift lj ~^ lj — Vs{r) ■ k, where Vs{r) = 2^$ is 
the superfluid velocity and P is the winding angle around 
a single vortex. As a result, contributions proportional 
to the magnetic field have been predicted to appear in 
thermodynamic and transport properties of line node 
superconductorsii^i^i^iii^i^iiSiiii. 

Magnetic contributions that appear because of 
the Doppler shift will strongly depend on the po- 
sition of the magnetic field with respect to the 
nodes. Consequently, oscillations in the DOS and 
related quantities have been predicted for super- 
conductors with line nodes in a rotating magnetic 
field,^'i^^'ii^'i^^'i^^^^^^ and have been ob- 
served experimentally in in-plane thermal conductivity 
in YBa2Cu307] ^^i^^ Similar results have been found for 
other unconventional superconductorS )^^'^^i^^'^^i'^°'^^i'^^ 
including YNi2B2G^ and PrOs4Sbi2,— which are re- 
ported to have point nodes instead of line nodes. Os- 
cillations in the angular field dependence of the specific 
heat have also been observed in several unconventional 
superconductorsi'^^i'^^'^^i'^^'^^i'^^'^^'^" 

Volovik's proof that, in superconductors with line 
nodes, delocalised quasi-particles have a greater contribu- 
tion to the low-energy DOS than vortex localised quasi- 
particles does not extend to superconductors with point 
nodes. However, we are interested in finding the oscil- 
latory component of the DOS in a rotating magnetic 
field, for which the semi-classical method, applied to de- 
localised quasi-particles, may be valid. For delocalised 
states, oscillations are obtained just by Doppler-shifting 



the quasiparticle energies. For localised states, the am- 
plitude of oscillations may be found by calculating sep- 
arately the DOS for the case when the field points in 
the direction of the nodes and the case when the field is 
perpendicular to the nodes and subtracting the results 
for the two cases. The former case corresponds to the s- 
wave result, which gives the contribution from localised 
quasi-particles as iVs_ioc. ~ Np^'^/B? ~ NfH/Hc2, 
where £, is the coherence length and R is the inter-vortex 
spacing. This sets a lower bound on the DOS contri- 
butions from localised quasi-particles in superconduc- 
tors with line nodes or point nodes. The upper bound 
of the DOS contribution from localised quasiparticles 
in superconductors with nodes is found when the field 
is parallel to the nodes. Volovik found that for su- 
perconductors with line nodes, the localised quasipar- 
ticles contribute Mine-ioc. ~ Nf^/hJh^/ log y/H\^jH 
which he found to be less than t he delocalised contri- 
bution Mine-deioc. ~ Np^J H/ Hc2- The osciUatory con- 
tribution for the delocalised states is contained within 
Minc-dcioc, whilc the oscillation amplitude for the lo- 
cahsed states is ^ iVune-ioc. - iVs-ioc. « -/Vuno-ioc- Thus 
the DOS oscillations in a line node superconductor are 
dominated by the delocalised contribution, and the semi- 
classical treatment is valid. For a point node super- 
conductor, one finds that the delocalised contribution 
is iVpoint-dcioc. ^ A'^point-ioc. ~ A^s-ioc- Again, the os- 
cillatory contribution for delocalised states is contained 
within A^point-dcioc. while the oscillation amplitude for 
the localised states is found by comparing A^point-ioc. to 
Ns-\oc.- This suggests that while both localised and de- 
localised states contribute to the DOS in point node su- 
perconductors, the oscillatory component is dominated 
by the delocalised states. We will assume that this is 
the case, but a more thorough investigation of the role 
of vortex localised quasiparticles in point node supercon- 
ductors is warranted. 

Field dependent thermal conductivity measurements 
are usually performed in one of two experimental con- 
figurations. In layered compounds in which the c-axis 
conductivity is low, such as in YBa2Cu307, the in-plane 
conductivity is usually measured with the B-field rotat- 
ing in the same plane. Then one in-plane component of 



the current will be parallel to vortices produced by the 
B-field while the other in-plane component of the current 
is perpendicular to the vortices. This introduces compli- 
cations when calculating the different components of the 
current averaged over the vortex lattice, since a different 
kind of averaging procedure should be used depending 
on whether the heat current is parallel or perpendicular 
to the vortices;^ As a consequence, for a field rotating 
in the xy plane, in-plane components of the conductiv- 
ity will oscillate with a period of twice the field angle 
even when there are no nodes at all, and this oscillation 
will dominate any nodal contributioni^i^ Thus, when- 
ever possible, the preferred set-up is to measure the heat 
currents perpendicular to the rotating B-field.— 



II. DENSITY OF STATES AND THERMAL 

CONDUCTIVITY FOR SUPERCONDUCTORS 

WITH LINE NODES 



In this section we consider a d-wave superconductor 
(such as YBa2Cu307) with line nodes along the direc- 



tions kx 



±ky and a magnetic field applied in the 



fcjjfcy-plane at an angle e with respect to the x axis. In 
the vicinity of a node, the gap function takes the form 
A(A;) « Vgk2, where fe2 points perpendicular to the node 



in the a;y-plane and Vg 



dAjk) 
dk 



is the gap velocity. 



The quasiparticle energy is E(k) = yje'^{k) -\- A2(fc) « 
v'^pk\ +Vgk2, where fci points in the direction of the 

node. Thus, in the vicinity of a node, the Green's func- 
tion takes the form 



In the superconductor PrOs4Sbi2, the 

pairing symmetry is widely thought to be 
unconventional,3"i4i,42i43^^^^^^ ^i^j^ gpi^ 

triplet pairing^ and broken time reversal symmetry.'*! 
Power law behavior has been observed in many 
thermodynamic and transport measurements at low 
temperature,— i^'^i^'^ which suggests the existence of 
nodes in the gap function; however a nodeless gap func- 
tion has been observed in some experiments i^°i^*'^^i^^ 
Oscillations of the thermal conductivity in a rotating 
magnetic field are another indication that there are 
nodes in the gap function.-- In previous works^^^*^ 
we have attempted to determine the symmetry of the 
superconducting state in PrOs4Sbi2 using available ex- 
perimental results. Among the various possible choices, 
we selected the spin triplet paired states belonging to 
the three dimensional irreducible representation T„ of 
the point group Th with symmetry D2{C2) x /C and 
order parameter components (0,0,1) and D2{E) with 
components (i|?72|, 0, |?7i|). We label these phases 'A' 
and 'B' respectively. The A phase is unitary, and has 
two cusp point nodes in the directions ±[0,0,1], while 
the B phase is nonunitary and has four cusp point nodes 
in the directions [±sin0OjO,± cos0o]i where (/)o is an 
angle determined from phenomenological parameters.— 
We will consider these and all other symmetry-allowed 
phases with point nodes. 



G(fc,iw„,r) 



[iCJn + aj(T"))^ + v'^pk\ + v'^kl 



(1) 



where iJ)„ = iujn + i^a, OLjif) = Vs{r)-kpj is the Doppler 
shift at the jth node and Fq = — 5Sret(w = 0) is the scat- 
tering rate at zero energy. The self-energy E is derived 
from the T-matrix formalism for impurity scattering and 
is the solution to the self-consistent equation^i^ 






(2) 



where F is proportional to the impurity concentration, c 
is related to the the phase shift (5o, c = cot (5o and 



G, 



o(iw„, r) = —— ^ G(fc, ia)„, r). 



nNp 



(3) 



Np is the density of states at the Fermi surface. In the 
unitary limit (c = 0) this yields a self-consistent equation 
for the scattering rate^i^ 



Fq = ir'^NpVFVgT 



In 



pI 



To 



ai{ry 



^{ai{r)+Tl){ai{r)+Tl) 



'Mr^,^^-,(Mr) 



V To 



■(4) 



In this article we calculate the DOS and residual trans- 
port under an applied magnetic field. We consider both 
the clean and dirty limits, in which the impurity scatter- 
ing rate is much smaller or greater than the Doppler shift, 
respectively. For the purpose of comparison, we begin by 
stating in Section II results for the residual DOS and 
thermal conductivity for the d-wave (line node) super- 
conductors. Section III is devoted to point nodes applied 
to PrOs4Sbi2. In Section IV we compare our results to 
experiment. Concluding remarks are made in Section V. 



where po is a cutoff and ai,2{r) are the Doppler shifts at 
two opposite nodes and can be written as 



ai(r) 
a2{r) 



2mr 
kp 

2mr 



sin (3 sin (7r/4 — e) (5) 

sin /3 cos (7r/4 — e) (6) 



where r is the distance from the centre of the vortex core, 
/3 is the vortex winding angle and e is the angle of the 
magnetic field relative to the x-axis. 



A. Density of States 



The DOS is given by 



N{uj,r) 



1 






^Gret{k,UJ,r) 



(7) 



where the integral over k is evaluated as the sum of four 
separate volume integrations centred about each nodei^ 
Then the DOS at the Fermi energy is 



iV(0,r) 



To 



TT VpVg 

ai{r) 



In 



pI 



Tn 



X tan 



tan' 



_i fai{r)\ ^ a2{r) 



Tn 






In the clean limit (ro/|a(r)|) — > and 
|ai(r)| + |a2(r)| 



N{0,r) 



2'KVpVg 



2mr 



sin/3| 



max [| sme|, | cose| 



2%/2' 



TTVFVg 



(8) 

(9) 
(10) 



This result necessarily has the same form as the finite 
frequency DOS N{uj) ~ \uj\ of superconductor with line 
nodes. Averaging over the vortex unit cell we obtain the 
result as in Ref. Il3, 



(^(0,r)>^ 



1 



7ri?2 



R r2-n 

drr I df3N{0,r) (11) 



?0 







AfFi^max[|sine|,|cose|] (12) 

R 



where ^ ~ \/w7- Evidently, there are four- fold oscilla- 
tions in the DOS as a function of the field angle e. 
In the dirty limit |a(r)|/ro ^ 1 wc find 



N{0,r) = 



To 



TT^VpVg 



2 In 



al{r)+al{r) 



^0 



(13) 



The first term is just the impurity induced DOS A^(0) 
so 



59 



SN{0,r)=N{0,r)~N{0) 



and the average DOS is 



1 



Vp 



An^r^ VgT() 



im0.r)),^Np^^^lJ^ 



sin^ 13 (14) 



(15) 



where Aq = t^ ~ NpvpVg. In this case, impurities 
remove the field directional dependence of the DOS. 
B. Thermal conductivity 



The thermal conductivity tensor is defined by the 
Kubo formulai^ In the limit T ^ it is expressed 
in terms of the imaginary part of the Green's function 

g^g59,60 

^^ - ^E^^^^ Tr[3G,et(0,r)3G,et(0,r)], 
^ -^ k 

(16) 
where ks is the Boltzmann constant and vp is the Fermi 
velocity in the direction of k. By again dividing the inte- 
gration over k into regions centred over each node, this 
eventually leads to 



^(0,r-) 
T 



kiE_,..vp.p r\,r,^^ 



3 {2TTYvFVg Jq 



2rg 



[{a,{r)+p)^ + rl] 



(17) 



where the integration variable isp — \/vpkf + f ?fcf . Per- 
forming the integration yields 



/t(0,r) ^ fc| Ej^i VFVF / ^ Qj(r) 

T 3 TTVpVg \ Fq 



tan 



-1 f(^jir) 



To J 2 
(18) 

where now vp is evaluated at each node. Summing over 
nodes yields 



jK(0,r) 



K(0,r)-K(0,0) _ k^vp 



T 

oi(t-) 

To 
ai(r) 



tan 



tan 




tan 



a2ir) 
To 



tan 



To y To y Fn y To 

In the clean limit a{r) 3> Fq, the thermal conductivity is 



— t;; — tan I T^ 



-1 I 0'2(r) \ ai(r) 



tan 




5k(0, r) _klvp ^ \^^{r)\ + \'^2{r)\ \^,{r)\~\c.,{r)\ 



T 



12 Wo 



To To 

|ai(r)|-|a2('')| |ai (r)| + |a2 (r-)] 

To To 



(19) 



(20) 



and the average over the vortex unit cell is 



T 



H 



2 VF Ao 




V2 



max sine, cose 



I sin (7r/4 - e) | - | cos (tt/A - e) | 



I sin (7r/4 - e) | - | cos (7r/4 - e) | 
I 



V2 



max sine , cose 



(21) 



In the dirty limit a(r) <C Fq we find 



SK{0,r) k'^vp 



"l(^')+"2(^) Ql('^)-"2('•) 
f5 f5 



fivr W„ 1 "l('')-"2(^) "l('')+"2('') 

9 \ p2 W2 

^ n ^ n 



The average over the vortex unit ceU is 



(22) 



/SKiO,r) 



H 



eJ-^^JLjHc.M 1 



VgTlH,2 \Hj\-sm2e 1 



sin 2e 



23) 



Thus, as in the DOS, impurities remove oscillations due 
to nodes in the diagonal components of the thermal con- 
ductivity. 



III. DENSITY OF STATES AND THERMAL 

CONDUCTIVITY FOR A POINT NODE 

SUPERCONDUCTOR 

We will first assume that there are an arbitrary number 
of linear (i.e. vanishing linearly with momentum) point 
nodes in the gap function, and that the gap velocities Vg 
are equal and isotropic around each node. We begin by 
finding a self-consistent equation for the scattering rate 
To analogous to Eq. 01 For point nodes, Eq. [3]is 



Go(0,r) = 



1 



E 



nodes 



2tt 



TtNf (2TT)^VFVg Jq 

—aj (r ) + iTo + p cos ( 
^ (-aj(r)+iFo)2-p2 



^6* sin 61 
Jo 



Po 



dp p^ 

(24) 



Vpkf + Vg(k2 + k^) and ki is parallel to the node 

while fc2,3 are perpendicular to the node. In Eq. 1241 we 
have again divided the volume of integration into parts 
each centred around a node. The integrations yield 



where the integration variable is 



Go(0,r) = 



Nf27T^Vfv'^ ^-^ 
" nodes 



-{To + iajir)f 



{To + iaj{r))pQ 



(25) 



Now we assume that there are four nodes which occur 
in pairs on opposite sides of the Fermi surface. Partners 
in each pair produce equal and opposite Doppler shifts. 



Summing over nodes we find 



Go(0,r) 



NF'2'^^VFVq 



4poTo+TT{al{r)+al{r)~2rl) 



T 

iVo 



(26) 



This result can easily be generalized to include more pairs 
of nodes. Equating the imaginary parts of (^5)1 yields the 
self-consistent equation for the scattering rate Fq, 



Fn = 



NfVfvJT 



2 p,ro-^rl + ^{al{r)+a^,{r)y 



(27) 



This equation describes how the scattering rate due to 
impurities is modified in the presence of Doppler shifted 
quasiparticles. 

As in Section II, we will assume that the magnetic field 
is parallel to the xy-plane with an angle e from the x axis. 



H = H{cos e X + sin e y) . 
The supercurrent is 
1 



(28) 



Vsir) 



2mr 



(— sin e cos /3x -I- cos e cos (3 y + sin (32) .(29) 



For now we will assume that all pairs of nodes are in the 
xy-plane at the positions 

fefi — ±fci?(cos(/>oa: — sin(/)oy) (30) 

kF2 = :izk F {cos (pax + sin (l)oy). (31) 

The angle 00 is zero in the A phase of PrOs4Sbi2 (and 
the gap function is doubly degenerate) and (/)o f^ in the 
B phase. This corresponds to the choice of the domain 
(1,0,0) of the A phase and the domain (|?7i|, z|?72|, 0) of 
the B phase. In each phase, two other domains are pos- 
sible and these will be discussed in the next section. 
The Doppler shifts are 



ai(r) 
a2{r) 



±Vs{r) ■ kFi 
±- cos/3[- 



= ± 



2mr 

±Vs{r) ■ kF2 
kp 



2mr 



sin 00 cos e — cos (po sin e](32) 
cos /3 [sin 00 cos e — cos 00 sine] (33) 



The following averages over the vortex unit cell will be 
useful: 



{alir)-al{r))^ 



4ln 
-l-sin^ 






— ) [cos^ 00 sin^ e 
Co, 



cos e 
R 

To 



sin 200 sin2e.(34) 



A. Density of states 




The DOS is given by Eq.[7l Using © and (|26|) we find 


r,p2 

N{Q,r) - -^^ 


-PO , -1 fPo\ 
tan — 

.To \^oJ 


2^0 \J^o/J 


■ (35) 



FIG. 1: (Color online) Oscillations in the density of states 
(Eq. I37p for different values of (j>o as a function of the field 
angle e. The dotted line is for (/jq = (A phase), the dashed 



PO 

line for (j>o = arcsin(0.3) and the bold line for 
7&ln(%^). 



N' 



Nf 



In zero magnetic field we retain our previous result for the 
impurity induced density of states^, then the magnetic 
contribution is 



B. Thermal Conductivity 



SN{0,r) 



al{r)+al{r) 



(36) 



which depends upon the Doppler shifts to the same power 
as that of the frequency in the low frequency DOS in 
zero magnetic field, N(ijj) ~ w^ for superconductors with 
point nodes. Taking the average over the vortex unit cell, 
we get 






[cos^ 4>o SIT? e+sin^ (po cos^ e] . where the 



Nf 

(37) 
Thus we find that the DOS oscillates with rotating mag- 
netic field as cos2e and is universal i.e. independent of 
the scattering rate. 



Beginning with Eq. 1161 we divide the volume of inte- 
gration into parts centred around each node, 



K(0,r) 
T 






2^1=1 



VpVp 



{2tt)^vfv^ 



PO 



dp p 



2Tr 



J- 



^6* sin 61 



0.8 







1.0 



integration variable is 



_r2i2 

oJ, 
again p 



(38) 



Vpkf + Vg{k2 + ^D- The integrations yield 



K(0,r) 
T 



fc| Y^j^l^FVp 



127r2 



VpVo 



To 



tan' 



-1 "iM Q'3(^)ro 



«(r) 



^0 



(39) 



where again our previously derived expression for the 
residual conductivity in zero magnetic freld^*^ is recov- 
ered. The matrix vpVp for one node is equal to the con- 
tribution for the node on the opposite side of the Fermi 
surface, but Uj (r) changes sign at opposite nodes, there- 
fore terms which are odd in aj{r) will vanish. The sum 
over nodes yields (keeping only the magnetic part) 



T 



Kg 



1 



127r vpv'gTQ 



/r2 



Vp 



12^ v^To 



[a\{r){vFVp)i + al{r){vpvp)2] 

{a\{r) + al{r)) coa^ (j)o i(a^(r) - af (r)) sin2(/)o O' 

i(a2(r)-a2(^))gin20o {al{r) + al{r)) siv? cj)^ 





(40) 



Finally we perforin the average over the vortex unit cell. 



SK{0,r) 
f 



H 



4^— In 




cos^ (j)Q [cos^ 00 sin^ e + sin (f>fj cos^ e] 

— I; sin2(/)o sin2e 





— J sin 2(j)o sin 2e 
sin^ (/)o [cos^ 00 sin^ e + sin^ 0o cos^ e] 




(41) 



where Aq 



Co 



NpVFVg. 



corresponds to 



The A phase of PrOs4Sbi2 
0, and the only component of the 



thermal conductivity which is non-vanishing is Kxx 



2 

sm e. 



Other domains 

The phase D2{E) has two other nodal configurations^^ 
which may be found by applying the operation C3 on the 
components (|f7i|, i|?72|, 0) or directly on the gap function. 
The second domain we consider is when the nodes are in 
the fc-cfc^-plane. Then the B phase has order parameter 
components (i|r72|,0, |r7i|) and the A phase has compo- 
nents (0, 0, 1). Then the positions of the nodes are 



kpi = ±fci?(— sin0oi + cos0o^) 
kp2 ~ ±kp(s\ii4>QX + cos 4>qz) 

and the Doppler shifts are 
kp 



(42) 
(43) 



ai{r) 
a2{r) 



±- 



± 



2mr 
kp 

2mr 



sin (f>o sin e cos f3 + cos 0o sin f3] (44) 
— sin 00 sin e cos /? -I- cos 0o sin /3](45) 



In the average over the vortex unit cell (aj(r) — a2(r)) 
vanishes, and 



{al{r) + al(r))^ 






[sm 00 sm e-|-cos 0oJ. 
(46) 



Then the B phase thermal conductivity is 



f 



e^^ln 



H 




• 2 / -2 

sm 00 sm e 



cos^ 00 
and the A phase thermal conductivity is Hzz ~ constant 

In t: 

plane 



In the third domain the nodes are found in the kykz- 



kpi 
kp2 



±fcF(cos0oj/ — sin0oz) (48) 

±fci?(cos0oy -f sin0oS) (49) 

In the average over the vortex unit cell, a\{r) -— a\{r) 
again vanishes, and we find 



i?2 



H, 



{al{r) + al{r))^ ^ _Z in ( _££ ) [sin^ 0o-fcos2 0o cos^ e] 

(50) 



H 




FIG. 2: (Color online) k^z (Eq. |47|| for different values of 0o 
as a function of the field angle e. The dotted line is for 0o = 
(A phase), the dashed line for 0o = arcsin(0.3) and the bold 

one is for 0o = f . (^L/T) ~ fc|^^ln(%^). 
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FIG. 3: (Color online) Kzz (Eq. I51|l for different values of 0o 
as a function of the field angle e. The dotted line is for 0o = 
(A phase), the dashed line for 0o = arcsin(0.3) and the bold 



(47) one is for 0o = f . {k'^z/T) ~ k 



2 ufA, 



■In I 



Bi;^ H^2 



Then the B phase thermal conductivity is 



/5K(0,r) 



T 



2VF^l H 



H 



V^To Hc2 



In 



Hc2 

IT 




cos^ 00 cos^ e] 



(51) 



and the only non- vanishing component in the A phase is 



f^yy ^^^ COS ^C 



Domain averaging 

In real situations, one may expect that either a single 
domain will form either because of sample shape or ap- 
plied strains or fields, or that all three domains will be 
present. If all three domains are present then detailed 
knowledge of the domain structure is required to calcu- 



late the conductivity. Lacking such knowledge, we con- 
sider two limiting cases: i) serial domains and ii) parallel 
domains. When the domains are in series the conduc- 
tivity is K = (k^ -|- k^ -|- k'^ ) which vanishes in all 
components. When the domains are in parallel the three 
conductivities are simply added: 



J 



f 



sin^ el 






I sin2 20, 



'0 



H 



4 sin' 200 




sin2e 



i sin^ 200 



cos^ e 



-isin 

;i-fsin2 20o) 




)o sin 2e 

i sin^ 200 




(52) 



r 



for the B phase, while the result for the A phase (0o = 0) 



SK{0,r) 
T 



H 



sin^ e 
cos^e 
1 



(53) 



Also, the domain averaged density of states is constant 
(has no oscillations). 



Other nodal configurations 



According to Table I of Ref . Im^ there are other nodal 
configurations corresponding to other superconducting 
phases which should be considered. Superconducting 
phases with cusp point nodes in tetrahedral supercon- 
ductors are summarised in Table [H 



Nodes 


Symmetry 


IR 


channel 


4 nodes [0, cos 0o, ± sin 0()] 


D2{E) 


r„ 


triplet 


2 nodes [0, 0, 1] 
6 nodes [0, 0, 1] 


D2{C2) X IC, C-i{E) 
Cs X IC, C^{E), /C, E 




triplet 
singlet 


2 nodes [1,1,1] 
8 nodes [1,1,1] 
8 nodes [1,1,1] 


CiiE) X /C 
T{D2), D2 X /C, D2 
T{D2) 


Eg 

Eu 


triplet 
singlet 
triplet 


6 nodes [0, 0, 1] and 
2 nodes [1,1,1] 


Cs{E) 


T, 


singlet 



TABLE I: Cusp point nodal configurations, their associated 
symmetries, order parameters labeled by irreducible represen- 
tation (IR) and pairing channels for a tetrahedral supercon- 
ductor (after Ref. [13). 

For eight point nodes in the [111] directions, the 

thermal conductivity is 



(5K(0,r) 
T 



H 



1 


2 


-sin2e 


- 


- sin 2e 


2 


\ 





2 



(54) 



The thermal conductivity for two point nodes in the 
[111] directions is 



5K(0,r) 
T 



H 



sin2e\ 


/111 
1 1 1 

V 1 1 1 


2 ) 



(55) 



Such a phase has three other domains, each with a single 
pair of nodes in the directions [1 — 1 — 1], [—11 — 1] and 
[—1 — 11]; the parallel domain averaged conductivity is 
equivalent to ([5i)) . 

The thermal conductivity for six point nodes in the 
[100] directions is given by ([55]) . 

The thermal conductivity for six point nodes in the 
[100] directions and two point nodes in the [111] 
directions is given by the sum of (|53p and ([SS]) . Such 
a phase has four domains; the parallel domain averaged 
conductivity is given by the sum of ((53| and (jM)) . 

Eqs. HniJ71 and [5T]|55l are summarized in Tablc|TTl In all 
cases, the highest harmonics which appear are two-fold 
oscillations, stemming from the fact that contributions 
from pairs of nodes are additive and proportional to the 
square of the Doppler shift. 



IV. DISCUSSION 

So far there has only been one report of thermal con- 
ductivity in a rotating magnetic field, namely the results 
by Izawa et al.r^ who measured Kzz and found four- fold 
oscillations near Hc2 and a sharp transition to two-fold 
oscillations at a lower field. We do not obtain four-fold 
oscillations for any of the point node configurations we 
considered and so we conclude that the formalism we 
have used is inapplicable in large magnetic fields. One 
possible source of error is that we have omitted contribu- 
tions from quasi-particle states localised in vortex cores, 
and that these states may dominate the oscillatory con- 
tribution to the density of states as the field increases and 
the vortices become closer together. Another possibility 
is higher order in a(r) (the Doppler shift) contributions 
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Nodes 


/^xx 


Kyy 


f^xy 


f^xz^yz 


Kzz 


4 nodes [cos <^o, ± sin 00,0] 






C 


c 


S 








4 nodes [± sin ^o , 0, cos (jio] 






C 











c 


4 nodes [0, cos (j>o, i sin (f>o\ 









c 








c 


domain average 






c 


c 


s 





1 


2 nodes [1,0,0] 






c 














2 nodes [0,1,0] 









c 











2 nodes [0, 0, 1] 


















1 


domain average/ 6 nodes [1, 


0, 


0] 


c 


c 








1 


2 nodes [1,1,1] 






s 


s 


s 


s 


s 


domain average/ 8 nodes [1, 


1, 


11 


1 


1 


s 





1 


6 nodes [1,0,0] and 2 nodes 


[1 


,1,1] 


c+s 


c+s 


s 


s 


s 


domain average 






c 


c 


s 





1 



TABLE II: Oscillatory contributions to the thermal conduc- 
tivity with a field rotating in the xy plane for various nodal 
configurations, 's' stands for sin2e, 'c' stands for cos2e, '1' 
stands for no oscillations and '0' means that the component 
vanishes, e is the angle of the field with respect to the x axis. 



become important as the field is increased. We do not 
obtain four-fold oscillations simply because we did not 
retain contributions to the density of states and thermal 
conductivity for powers of a{r) higher than two. In any 
case, unlike dx2-y2 line node superconductors, the four- 
fold oscillations reported in Ref. 3u are not related in any 
simple way to the nodal structure of PrOs4Sbi2. 

Our results may be applicable to the lower field ther- 
mal conductivity measurements in which two-fold oscil- 
lations are found. The lower inset of Fig. 1 of Ref. 'SO 
shows a nearly linear dependence of Kzz on H, in rough 
agreement with HlogH dependence expected for point 
nodes (see Eq.HT]). Fig. 2b) of Ref. l30l clearly shows two- 
fold oscillations of the form k^z ~ cos 2e and not sin 2e. 
This indicates that the most likely superconducting phase 
of PrOs4Sbi2 is D2{E) which belongs to the three di- 
mensional order parameter T„ and that a single do- 
main with order parameter components (0, It^iJ, i|7y2|) or 
{i\r]2\,0, \r]i\) was measured in Ref. ISOl We note that this 
phase agrees with various properties observed in other 
experiments, including triplet pairing,'*^ broken time re- 
versal symmetry'*^ and broken C3 symmetry4^i^ 

Sakakibara et al.''*' measured the field-angle dependent 
specific heat and found four- fold oscillations, but unlike 
Izawa et al., no four- fold to two- fold transition was ob- 
served. In this case we must also attribute the four- 
fold oscillations to corrections beyond the semi-classical 
methods we have used. Sakakibara et al4S verified that, 
in their set-up, superconductivity has no preferred orien- 
tation, implying that their results are domain averaged. 
Since two-fold oscillations are not observed at all, these 
results may be consistent with any of the domain aver- 
aged configurations shown in Table |TT1 

Similarly to computational issues pertaining to domain 
averaging, the thermal conductivity must also be aver- 



aged over the vortex lattice. In all of our calculations 
we performed the vortex average as a simple areal av- 
erage over a plane perpendicular to a vortex, as shown 
in Eq. 1111 This procedure is appropriate when the heat 
current is parallel to the vortices. For currents in other 
directions a different averaging procedure should be used, 
which results in a more complicated field dependence 
of the oscillation amplitudes than what we have shown 
here. The correct procedure is an average of n over 
paths through the vortex lattice, which is in fact more 
involved than the series average (k^^)^^ described in 
Ref. 7. This means that the vortex averaging calcula- 
tion for any in-plane component of the conductivity will 
vary with the field angle, producing oscillations ^ cos 2e 
which are unrelated to nodes and which will dominate 
over any nodal contributions. ^^•^'^ Thus observations of 
oscillations '^ cos2e in k^x, Hyy or in off-diagonal com- 
ponents of K measured with an in-plane current should 
not be interpreted as evidence of nodes. 



For Kz 



and Kyz, when the current is perpendicular 



to vortices, there will not be oscillations due to the vor- 
tex averaging. If observed, oscillations may therefore be 
attributed to nodes. A small oscillatory contribution will 
arise from mixing with the other components of n via the 
vortex averaging procedure (which in general does involve 
averaging k~^) but we expect it to be small compared to 
the oscillations originating from nodes. 



V. SUMMARY 



We have reviewed previous works concerning field- 
angle dependent DOS and thermal conductivity for line 
node superconductors using a semi-classical method, and 
applied the same method to point node superconductors. 
This method neglects vortex localised quasi-particles and 
retains only the contribution from extended, nodal quasi- 
particles to the density of states. Clearly there are limita- 
tions to this approach; in particular it cannot be expected 
to produce an accurate estimate of the total low-energy 
density of states in point node superconductors. However 
it may be a reasonable way to estimate the field-angle 
dependent oscillatory component of the density of states 
and related quantities for fields Hd < iJ <C i?c2- We 
find that in point node superconductors there is no dif- 
ference between the clean and dirty limits, unlike in line 
node superconductors in which the different limits pro- 
duce significantly different expressions for the oscillatory 
part of the thermal conductivity. Considering all possible 
configurations of point nodes in a tetrahedral supercon- 
ductor, we find that the superconducting phase D2{E), 
which we previously proposed based on other experimen- 
tal evidence, best accounts for field-angle dependent os- 
cillations in the thermal conductivity of PrOs4Sbi2. 
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